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Summary 

A formalism is described for micromechanic impact of composites. The formalism consists of numer- 
ous equations which describe all aspects of impact from impactor and composite conditions to impact 
contact, damage progression, and penetration or containment. The formalism is based on through-the- 
thickness displacement increments simulation which makes it convenient to track local damage in terms 
of microfailure modes and their respective characteristics. A flow chart is provided to cast the formalism 
(numerous equations) into a computer code for embedment in composite mechanic codes and/or finite ele- 
ment composite structural analysis. 


Introduction 

Composites in aerospace vehicles in general and composites in engine structures in particular are 
subjected to inadvertent and design-for impact loads. References [1-3] describe various types of impacts, 
approximate methods, while refs. [4-6] describe different experimental techniques to evaluate compos- 
ite impact. References [7,8] describe impact damage in composite structures by classical formulations 
while refs. [9,10] describe impact damage in composite engine structures. Damage tolerance and dynamic 
loading in composites are extensively described in two recent books [11,12], From those studies as well 
as others and continuing research on three-dimensional finite element simulations, it is fair to note that 
impact damage in composites is very complex and is mostly evaluated by testing [11-13]. 

The complexity arises from the multitude of composite fracture modes that are caused during impact 
and the attendant failure mechanisms associated with those fracture modes. Impact damage for specific 
designs can be evaluated by finite element simulations. However, local tracking of damage during im- 
pact and identification of local fracture modes is a very tedious and often an elusive process. Because the 
finite element model for composite impact is usually very large, it requires substantial personnel expertise 
and effort, and long computer times. Therefore it is cost prohibitive if not impossible to use finite element 
simulations in preliminary designs of composite structures for impact resistance. It is also fair to note that 
an inclusive micromechanics treatise of composite impact has not been developed or even proposed as 
yet [14-18]. The work reported in [14-15] is a step in that direction. 

It is timely then to develop a micromechanics based generic and inclusive treatise of composite im- 
pact. The objective of this report is to describe the development of the formalism for a micromechanics 
treatise of composite impact. The treatise covers all aspects of impact damage from contact to penetra- 
tion and includes cratering and spallation. It describes local damage as it occurs and the various fracture 
modes that are induced during the damage process. Strain rate effects are neglected because of the in situ 
three-dimensional constraints in the matrix. Only point impact is considered since it is assumed that this 
will be coupled with structural analysis for global response. Equations are derived for penetration time, 
residual velocity, rebound velocity and backface splitting. Equations are also derived for impact force- 
time trace and impact containment. A flow chart is outlined which cast the formalism into a computer 
module computational algorithm which can be made part of available composite mechanics codes and 
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general purpose finite element for composite structural analysis. Because of its inclusiveness and 
generality the treatise is approximate relative to mathematical rigor but consistent in the approximations 
with mechanics of materials approach. 


Nomenclature 


A 

D 

E 

F 

G 

8 

h 

K 

M 

R 


r 

S 

T 


t 


V: 


ind 


V 


x,y,z 

a 

P 

AT 

e 


0 

1 

v 

P 

2 
a 

{},[] 


area 

projector diameter 

energy, modulus, E.E. elastic energy 

force 

shear modulus 
gravitation constant 
crater height 
K.E. kinetic energy 
moisture, mass 

projectile radius, transformation array 

contact area radius 

strength 

time, tension, temperature 
thickness, time, t e ply thickness 
volume of indentation 
velocity 

projectile penetration depth 

coordinate axis 

thermal expansion coefficient 

moisture expansion coefficient 

temperature change 

limit strain 

strain 

angle contact area periphery (figure 1), ply orientation angle (figure 6) 

laminate property 

Poisson's ratio 

material density 

summation 

stress 

vector, array, respectively 


Subscripts 

C compression 

c laminate, contact crushing 

crz crz composite, rz direction 

F bending (flexural) fiber 

i indentation 

o initial property 

p projectile related, penetration 

S strength 

s shear 

x,y,z coordinate direction (figure 13) 

1,2,3 material axes (1-along the fiber) 
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Some Geometric Relationships 


Express the radius, r, and the angle, 0, of the contact area as a function of z the local indentation of an 

s R p (figure 1). 

(R„-zf + (rf=R 2 p 


impactor with a spherical tip of radius R p (figure 1). 


R 2 p -2R p z+z 2 +r 2 =R 2 p 


( 2 . 1 ) 


r 2 =2 R v z-Z 2 =z(2R p -z ) 


T / \Tt/2 

=[z(2R p -z)\ 


0 = sin 


f \ 
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= sin 1 • 


R n 

Rn 
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P 


(2.1a) 


( 2 . 2 ) 


The volume of the spherical cap with radius r and depth z is derived as follows: 
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(2.3) 


(2.4) 
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Penetration Description 


Penetration is a time dependent process. Here we consider very fast penetration that is the structure 
responds only locally by material damage (stress waves do not reach the boundaries). 

The material getting indented is equal to the volume of the spherical cap in the indentation (figure 2). 
The energy required to crush this volume of material is given by 

E c = \VindS c 33C (3-D 

Where E c denotes the energy of indentation, V jnd is the crushed volume, and S c33c is the flat-wise lami- 
nate strength. If we further assume that this equals the kinetic energy absorbed (relatively stiff indentor 
compared to composite flat-wise stiffness) a projectile with mass M and an initial velocity v 0 , then for 
energy conservations the kinetic energy of the projectile equals the energy expanded E c to crush the 
material: 


E= \ — M„v. 


P o 


2 ^ind R c33C 


Mp v o ^ind R c33C 


f Vi/2 

VindS c 33C 


M r 


It is noted that v is the projectile velocity prior to contacting the target. Check units: 


in? iff in? 

1/2 

in? if in? 

in? if)/ in? in. sec 2 


ijyj [in/ sec 2 j 


/ 9 ,xl/2 

v 0 = I in. /sec~J = in. /sec 


(3.2) 


Therefore, the units check. 

The next step is to assume that the projectile has first to crush the impacted surface up to its 
equator and beyond until the crushing force exceeds the remaining thickness shear-stress-resistance to 
push it through as a plug. 


v ind = ^ D p wh en (z^ R P ) ( 3 - 3 ) 

v i n d = -^ + n R p(z- R p) when (z>/? p ) ( 3 - 4 ) 

Equation (3.3) is half the spherical volume for a spherical projectile and equation (3.4) is that plus the 
volume of the shear plug figure 3. Then the crushing velocity required at z > R p and at z < R is given 
respectively by 
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±nDl + nRj(z-R p )]s c 


1 3 

-nDpPp/g 




2gz 2 (3R p -z)s c33c 


* 4 / P, 




Additional energy will be required to push the material ahead as a shear plug for penetration. That energy 
(E ) is given by assuming the friction against the projectile is part of the shear resistance, figure 3 and the 
projectile velocity is v and the time t s required during the shearing plug: 


“ 9 S C rzs[ nD P t s ) 


E s = - 
s 2 


r p* z p 


— M p v s 2 = — 
2 1 2 


9 ^crzS ( KD ph ) 


nDpt s S crz s 


Another way to think is v s is the projectile velocity the instant the shear block starts moving. For a 
spherical projectile of diameter D p and for M p > than the mass of the plug: 


n D p 

M p = -^( P /g) 


_ S^crzS 


rfPp 


il/2 f . 2 / - / 2\ ii /• 3 " 

m \in / sec Jib/ in 


in 1 {lb]ur' 


= — O.K. 


sec" I sec 


the units check. 

The total velocity required for penetration E p is calculated from the sum of the two energies 

F = F + F 

‘-‘p ‘-‘c T 


and the residual velocity is 


? ( 3z R p) S c33C | 3gt^ S crzS 
4R p Pp IRpPp 


It is noted that v p is velocity that the projectile must have to crush and then staid pushing the plug out. 
The residual velocity at the instant of v p is E r which is given by: 
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Ef — E a Ep — E a E c E s 


Dividing by M p the 
equation for the 
residual velocity is 


E r — M ,,v~ M p v c S crz ^nD p t s 


l p y o 


-il/2 




r 

-j 

> 

v r = 

2 2 r\ +2 

v o ~ v c — KDpt s 

ScrzS 

M n 



- 

P 


(3.9) 


v r is the residual velocity that drives the shear plug through, where v g is the projectile initial velocity, v c is 
the crushing velocity and the last term in the equation is the shear velocity. 

Several points are worth observing: 


'• v o=K' +y s] 

2. V n > l'’ 2 r + v2 J 


2 t 1/2 


2 1 1/2 


2 1 1/2 


v, = 0 


1/2 


V r = [v 2 0 - ( v 2 c + V 2 S )] 

Only partial penetration if at all. 


3. v 0 < |v z ( . + v z s ] 

4. For any velocity v Q and a given laminate with known S c3 3C and S crzS , a composite thickness t ( . 
can be found for projectile containment. 

5. Equivalently for a given t ( . and S c ^ c , an S a .„ s can be selected for containment. 

The equivalent force at the instant which the plug begins to move will be equal to the normal stress 
times the projected indentation surface area. For a spherical projectile at z = R p the surface area is: 


A :s — ' 


nD, 


and the force 


F c -A:sSc33C (3.10) 

F c =UD p 2 S c33c (3.11) 


This force or ratio thereof can then be used to perform evaluations for delaminations, splits and other 
fracture modes in the remaining laminate. 

It is important to note that S crzS very likely includes contributions from fiber fractures either in shear 
or in local bending. This implies that penetration will be dependent on the laminate configuration. It is 
suspected that delaminations and/or fiber crushing will occur when: 


1 2r 

A- — o c 3 3 C ^ n D p zS crz s 

_ i 2 1 7 

F c . - - TlDp S c 33c > -ScrFtf, 


[ interply strength! 
[at that inteiply j 
[Ply bending strength! 
[at that intraply j 


(3.11a) 


If fibers break in local bending the crushing force is 
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Where S crF is the flexural strength and t, is the thickness of that ply, respectively. 


Time to Penetration 

Time to penetration (T p ) is the time required for the projectile to penetrate the laminate. That time 
can be estimated by assuming that the following relationship holds in an average sense: 

1 T , 

2 v P T P =t P 

2t v 

T p = — (3.12) 

where v p is given in eq. (3.8) and t is the penetration thickness (t = z). Since both v Q and v p can be 
measured in a high-velocity-impact test, suitable tests can be conducted to verify eqs. (3.8) and (3.12) and 
indirectly eqs. (3.10), to (3.11b). It is noted here that equation (3.12) is only a rough estimate. 


Cratering 

Cratering occurs at the periphery of the contact as the material is pushed out due to high in-plane 
compressive stresses which are induced by the projectile ingress during crushing. The cratering is esti- 
mated by the following relationship (figure 4): 


£ 


f = —V £ 

°czz v crz err 


= —V 

czz v crz 


A33C 

E c33 


(3.13) 


The crater height is determined from 


crater ^ czz 


(3.14) 


where t , v crz , S c 33C , and e c33 are evaluated at the depth of crushing at that time. Implicit in eq. (3.14) is 
that no cratering will occur once shearing started. The corresponding time is approximated: 


At-, 

2 vJc tp 

2 f 

T y_ 


(3.15) 
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where v c is given by eq. (3.5). We can also determine whether spallation will occur as a result of cratering 
since we know the strain e ( .„ at any (t - z). 

A condition for estimating spallation is given by 

a c33 > S c33T 


where s c33T is given by 


'-c33^c33 




'c33C ' 


\ 33 
-c33 


= v S 
v crz ■ 


c33C > S c33T 


(3.16) 


Delamination and Back Surface Splitting 

These two conditions will occur if the penetration time is sufficiently long for shear waves to propa- 
gate in a region beyond projectile diameter. The shear stress waves velocities (v crz ) are given by 


f r ^ /2 
'fax. 


V Pc J 


(3.17) 


v cxy 


xl/2 


J cxy 


V l J c / 


(3.18) 


where G cn and p c are the corresponding laminate shear modulus and density, respectively. The distance 
these waves will travel at T c are estimated by 


? v crz^c 
X = y = v cxy T c > 
T c = 2 t c / v c 


(3.19) 


Recall that t is the crushed thickness and v the corresponding velocity. Using the last equation in the 

P C 

above two equations and also eqs. (3.17) and (3.18) and by assuming average wall stress 


where v p 


\l/2 
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J crz 


V He J 


2 t r 


V v c J 
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g(D 3 p + 3R 2 p (z-R p )) 
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-ill 2 


2Dlp 
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. 6t s ScrzS 
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(3.20) 
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Delamination will occur if: 


r > D p / 2 


an d ~D p S c33C > t s S crzS 


and back surface splitting will occur if in addition to eq. (3.21) the following conditions hold: 

d ran verse ' 

1 2 1 2 

- KD p S c33c >~h S c 22 T ply 

^splitting 


1 2 


longitudinal ' 


— nD p S c33 c >- t ‘f ^ c UF\ fiber 


breakage 


Note that the r in eq. (3.21) indicates an average radial stress wave velocity propagation. Values for r x and 
r are calculated from 


r v = y = 


T c cos 0 


T c sin 0 


where the x and y are the in-plane coordinates. These equations will give an elliptical wave propagation 
front, in general. 


Thickness Growth 

This growth is due to Poisson's effect. The ply in-plane unrestrained strains along the material axes 


Those in the in situ laminate are 


{e — {Vr/m } AT|cq> } +M{(3 m } 


{£i} = [%e]{e c } 


The ply in-plane restrained strains are the difference between eqs. (4.2) and (4.1) 


{ef} = [%e]{e c }-{e4 
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The thickness growth is obtained by multiplying eq. (4.3) by the appropriate ply Poisson's ratios 


- izz 



j = 1 and 2 


(4.4) 


The composite or laminate thickness strain then is obtained from 


j 

£ czz ~ ~ ^ {hfiizz ) 

i r . , 


c i = 1 


(4.5) 


and the laminate change in thickness is obtained from 


N f 

Af C =X( f « £ fez) ( 4 - 6 ) 

i = 1 

It is noted that these strains are in addition to those produced by impact. It should also be noted the 
ply in-plane strains are already known from laminate analysis. To obtain the thickness strains and thick- 
ness changes use the ply strains from laminate analysis in eq. (4.4) and proceed with eqs. (4.5) and (4.6). 


Elastic Indentation Stresses 

For an elastic indentation of depth Z, we assume continuous contact between projectile surface and 
target indentation. For this case, the in-plane and thickness-size local strains are uniform and the stresses 
will be proportional the local stiffness. Since the strains are uniform and the deformation conforms to the 
projectile surface, then at a depth Z it follows that the in-plane (e CTT and e^q) and thickness-wise (e ( „) 
strains are: 


_ A r 

£ err ~~ 

r 



£ 


czz 



(5.0) 


Note eq. (5.0) denotes completely recoverable indentation. Note also that for a spherical indentor at any z, 
r = \z.(2R jt - z)] 1/2 , eq. (2.1a) and q = sin -1 r/R p . Then, the in-plane ply strains along the material axis are 
given by (see figure 5) 




(5.1a) 
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£ 11 
e 22 


(5.1b) 


cos 2 0 sin 2 0 — sin 20 
2 

sin 2 0 cos 2 0 —sin 20 
2 

-sin 20 sin 20 cos 20 

The corresponding stresses are given by 

M = [£<#?} (5.1c) 

where E, is defined by: 

E m ~ v m E e.22 0 

l~ v £12 v £21 1 — V ^12 V T21 \Z(\\ 

vm E m E 122 0 • £ «2 ’ (5-ld) 

1 — v £12 v £21 1 — V ^12 V T21 l £ ^12. 

0 0 G tt 2 

Concisely the in-plane strains are given by 

{ £ c4=[%eU £ cr} (5-2) 

and the flat-wise 

{ e zf} = K/z] f Kz} ( 5 - 3 ) 

A more direct approach will be to obtain the three-dimensional strains from the three-dimensional 
transformation (see appendix A). 

{ e z^j = [^£/r,e]^{ e cz,z} (5-4) 

The last equation will provide all 6-strains from which the corresponding stresses can be obtained 
from the three-dimensional stress/strain relationships: 

{°zf} = [ £ c J{ e zf} (5-5) 


s = 

Equations (5.4) and (5.5) can be evaluated by incrementing eq. (5.3) by a single fiber cell 
and continuously checking fracture modes. Here s is the edge of a fiber square array, Ay is the fiber volume 
ration and is the fiber diameter. The strains fracture modes are obtained from 
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(5.6) 


and the stresses fracture modes from 


(5-7) 

If eqs. (5.6) or (5.7) are violated failure will occur by that fracture mode and permanent damage will be 
induced in the composite target at that location. Crushing will proceed as described in Penetration De- 
scription section. 

The elastic energy stored down to depth z will be given by 

E - E =\ X ( 5 - 8 ) 

/ Az=0 

where Az is the thickness increment and r is evaluated at z from eq. (2. la), in the section titled Some 
Geometric Relationship. 

The loss in kinetic energy then is 


■E-) Ioss =E.E. (5,9) 

when 

(5-10) 

No damage will occur and the projectile will bounce (rebound) with velocity equal to its initial velocity 
assuming no losses due to local heat or vibrations. To assure no heat and no vibration the contact time will 
have to be relatively small so that the change in local temperature is negligible: 


AT ~ 0 


(5.11) 


where the impact surface temperature and the shear stress waves have not propagated beyond the a 
distance t p (crushed thickness). The contact time for no damage can be evaluated from the rebound 
condition 


o 

1 

o 

II 

(5.12) 

2 z 


E c ~ 

(5.13) 

v o 



Assuming an initial velocity of 30.5 m/sec (100 ft/sec) and an elastic indentation of z ~ 0.013 cm 
(0.005 in.) without damage in any fracture mode 


T _ 2(0.005) 
c 1200 


0.010 

1200 


— x 10“ 5 sec = 8.33q.v 
1.2 
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a very small time indeed. 

Also the contact time for onset of damage can be calculated from eq. (5.8) by setting 

z 

E - E -(o ze >S ze ) = XL £ zf J{ a zf) dz ( 5 - 14 ) 

Az 

These simplified equations can be readily checked by using available experimental data [3-5, 12 and 
13]. Note that the aforementioned derivations are for hard object projectile. They need to be modified for 
soft impact. In that case the z increment will be on the projectile and r on the projectile flattening. 

The other important point to note is that the energy for elastic deformation must be added to that for 
crushing and penetration if it is rather significant, say, more than 3 percent of their sum. 


Force-Time Trace 

The force time trace can be determined for both elastic impact and penetration impact. Elastic Im- 
pact-elastic indentation the force at any local indentation z is obtained from (figure 6). 

F z = nr 2 o lzz (6.1) 


where r is given by 


r it/ 2 

r = \z(2 R p -z ) 


and o czz( is given by 


o 


tzz -~ E t 33 
l c 


and the corresponding time is given by 


T = 


2z 


V a +V z 


( 6 . 2 ) 


(6.3) 


(6.4) 


where v o is the projectile velocity at initial contact and v is determined from (K.E.—E.E.) is satisfied or 
equivalently from the equation: 


{m p v 2 z = 1 -M p vI~E. e .(z) 


(6.5) 


Unloading will be obtained by deincrementing z and the contact time will be obtained from z max 
which occurs when v = 0 or from eq. (6.5): 


E - E (z ma j) = ^ M p v l- 
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Impact Damage Response 


Impact damage response is evaluated by accounting for the energy expanded as impact damage 
progresses and accumulates (figure 7). The impact force is 


=nr 2 a czz 


(7.1) 


17 ( 2 ) _ _ 2 
t z ~ m ^c33C 


?(3) _ 


2jt Rr,t r S, 


p l c.crz.S 


(7.2) 

(7.3) 


d 4 ’= 0 


(7.4) 


The corresponding times are given by: 


T z l) =2z/{v 0 +v z ) 


(7.5) 


where — M n vZ = 

O F 


~M p Vg -E.E. 


(eq. (5.8), until C7 z€ =S €33C ); 


(7.6) 


T, (2) =2z/(v 0 +v,) 


(7.7) 


where — M ,.v? = 
2 1 


\m p vI-E 2 


(E 2 , eq. (3.2), 3 shearing takes place); 


(7.8) 


from the moment that crushing begins until the moment of transition from crushing to shear plug move- 

17(3) < f (2) 

ment, or from the moment that o 22 l ~ s f33 t0 ^ 1C moment z ~ z ; where z c is the penetration depth at 
the moment crushing ends; and v z is the projectile's velocity at that same moment. 


<3) 


= 2z c /(v 0 +v z )+ 2 ^ 

' > 1; 4-17 


V + V 
v z v z 


(7.9) 


where — M„vl = 
2 1 


\m p v~-E s 


(E s 3 due to shear) ; 


(7.10) 


(4) _ 2z c _ J _ 2 ( t c +R p Zc) 


■ 


- + - 


V a + v z 


v z +v e 


(7.11) 


where t is the composite target thickness and where v is the exit velocity obtained from equation (7.10). 
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Computation Simulation 


The penetration process can be simulated as described below. We can use the previous simple equa- 
tions to computationally simulate the penetration process: The block diagram in figure 8 is one way to 
implement the impact process described, by the equations in the previous sections. 


Typical Results 

Typical results were obtained by using the computer module generated from the block diagram 
(flowchart) in figure 8 in conjunction with a composite mechanics computer-programmed [19] and 
reprogrammed, respectively in JAVA computer language code. These results are from assuming a hybrid 
[0/90/0] s laminate where the 0°-plies are AS graphite-fiber epoxy-matrix and the 90°-intraply hybrid plies 
are AS graphite-fiber epoxy-matrix// S -glass-fibers/epoxy-matrix composites. The constituent properties 
used are typical for these composites. The results are for three different initial velocities 36.6 m/sec 
(120 ft/sec), 366 m/sec (1200 ft/sec), and 3660 m/sec (12000 ft/sec). The projectile is a 0.64 cm (0.25 in.) 
diameter steel ball. 

The force-time traces developed during impact and predicted by eqs. (7.1) through (7.4) are 
shown in figure 9. Three points are well worth observing: (1) the impact duration for each velocity, 

(2) projectile containment, and (3) laminate penetration. The locations for the projectile progression 
predicted by eq. (7.5) through (7.11) are shown in figure 10. Note containment for the two low velocities 
and early penetration for the highest. The velocity trace predicted by eq. (7.9) is shown in figure 11. Note 
the containment velocity is zero for the two low velocity impacts and very high residual (escape) velocity 
(about 90 percent retention) for the high velocity. The interply and intraply stresses develop during impact 
predicted by eqs. (3.1 la) and (6.3) are shown in figure 12. Note that the lowest velocity reached only 
69 Mpa (10 ksi) while both the others reached about 1 10 Mpa (16 ksi) denoting most likely interply and 
intraply delaminations for both of these. It is worth noting that the discontinuities in stresses shown in fig- 
ure 12 are from the assumption that the in-plane strains are continuous prior to damage. Suffice it to say 
that this type of the laminate local detail can only be obtained by a micromechanics impact treatise such 
as described herein. 


Conclusions 

An investigation was performed to develop an approximate formalism for describing composite 
damage due to impact. The formalism consists of numerous composite micromechanic equations describ- 
ing all aspects of impact from initial conditions (projectile, target) to containment and/or penetration. 
These equations describe damage initiation, damage progression and accumulation in terms of microfrac- 
ture and respective fracture modes. Equations are also included for impact-force/time history and time to 
penetration as well as conditions for back surface splitting. A flow chart is outlined to cast the formalism 
in a computer code module for embedment in composite mechanics codes. The formalism in essence 
constitutes an approximate treatise of micromechanics for composite impact. Results obtained from the 
computer module of this flow chart illustrate laminate local detail information that can be generated. The 
three velocities investigated may be characterized as follows: crushing initiation, extensive crushing, and 
hyper velocity impact. 
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Appendix A 


The three-dimensional strain transformations are given by the equations below where t, m, n are 
direction defined by the table below (figure Al). 


2 2 2 

£ cxx = h £ crr + m \ £ c00 + n \ £ czz 
2 2 2 

e cyy = h £ crr + m 2 £ cQQ +n 2 £ czz 

e CZZ — i\ £ 

err "f m 3 £ c00 n 2 £ czz 

£ cxy = 2t\£ 2 £ crr + ^ m \ m 2 £ cQQ + ^ n \ n 2 £ czz 
£ cyz = ^2^3 £ crr +2m 2 m 3 £ cQQ +2n 2 n 3 £ czz 
£ czx = 2 hh £ crr + 2m 3 m l e c00 + 2 «3 n l £ czz 


(Al) 


r, q, and z are attached to the projectile tangential plane and x, y, z are the laminate global axis. The 
undamaged ply stresses in the laminate at that point are obtained by: first, calculating the ply strains 
e tz7 , and e n -, along the material axis and second, by using the three-dimensional stress/strain relation- 
ships. The flat-wise (normal) ply stress is obtained by multiplying e by and the flat-wise ply 
shear stress from 


are obtained from 


i £ czz^B3)j 

° tyz ~ £ cyfilyz 
0 izx ~ £ ezxfiizx 


(A2) 


(A3) 


The interlaminar stresses are obtained from (assuming isotropic interply layer behavior) 



(A4) 


In eq. (A4) i denotes interlaminar and j denotes the jth interlaminar layer. 

Equations (Al) to (A4) are then compared to their corresponding strengths to identify local fractures 
and respective fracture modes. It is noted that eqs. (Al) to (A4) are general and can treat any material 
including isotropic hybrids and combinations. Thermal, moisture and residual stresses can be included 
as is done in the conventional laminate analysis. For local damage and progressive fracture the projectile 
induced strains are obtained by the incremental growth during the current time increment. 
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Figure 2. — Indentation diagram. 
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Figure 3. — Penetration diagram. 



Crater 
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Figure 6. — Impact force history diagram. 



Figure 7. — Impact damage response. 
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1) Input: 

Targe F.E. Model 
Impactor F.E. Model 
and initial cond 

T 

2) Select impact 
region and 
model as shown 

i J 

3) Use composite 
Mech to synthesize 
laminate properties 

E 

4) Initialize at contact 
T = 0, z = 0 

Ez = ra |x2/2 
h = 0 


5) Calculate 

dz = min {tc, Rp}/ 
n increments 

6) Increment z 

I 

7) Compute crushing 
force Fz 3.1 1 

I 

8) Integrate to get 
remaining energy 
Ez = Ez'-Fz dz 

~ 

9) If Ez < 0, set 
Ez = 0 

1 

10) Calculate 
Vind (Eg- 2.4) 

V 

11) Velocity: 

v = sqrt [2 Ez/Mp] 

V 

12) dT = dz/v 

I ~ 

13) Integrate time: 

T = T' + dt 

I 

14) Integrate crater 
height: h crater = 
h crater 1 + epslzz 
dz 

~ 

15) Test for failure 
modes 


16) If v = > 0 
and 

z < Rp loop 


Completed initial crushing in 
shape of projectile. Continue 
crushing with cylindrical walls. 



27) Test for failure 
modes 




30) Ez = Ez'-EzDrop 
Where EzDrop = 


Fz (thickness + Rp- 
z chg) (? see 3.6) 


If v is zero 


skip to 

print at end or 
z = thickness + Rp 


31) If Ez<0, Ez = 0, 
If v = 0, skip to 
print 


Figure 8. — Computational simulation flow chart. 
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0.00 1.00 2.00 3.00 4.00x1 0" 5 

Time, sec 

Figure 9. — Impact force developed in a [0/90/0] s 
interply hybrid laminate (90°-plies intraply hybrid 
graphite/S-glass). (1 lb = 4.45 N, 1 ft = 0.305 m). 



Figure 10 - Projective location during impact [0/90/0] s 
interply hybrid laminate. (90°-plies intraply hybrid 
graphite/S-glass). (1 in = 2.5 cm; 1 ft = 0.305 m). 
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I.fixio 4 



0.00 1.00 2.00 3.00 4.00X1 0" 5 

Time, sec 

Figure 1 1 . — Projectile velocity in a [0/90/0] s interply 
hybrid laminate during impact. (90°-plies intraply 
hybrid graphite/S-glass). (1 ft = 0.305 m). 


1.6X10 4 



Location, z, ft 

Figure 12. — Interlaminar and intralaminar shear 
stress developed in [0/90/0] s interply hybrid 
laminate during impact (90°-plies intraply hybrid 
graphite/S-glass). (1 psi = 6.9 Pa). 
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Figure A1 . — Coordinate transformation 
diagram. 
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